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INTRODUCTION 


Breaking of window glass due to heat from fires is a very commonly observed phenomenon. This 
phenomenon is of great importance as the glass breakage can be a cause of fire spread (if the window 
is between two adjoining rooms) or a broken window can act as a vent for the escape of toxic fire 
gases and an inlet for fresh air (if the window opens to the outside). Glass absorbs thermal radiation 
from the fire and also gets heated by convection from the hot air surrounding the fire leading to ther- 
mal stresses. The thermal stresses arise due to the temperature difference between the central portion of 
the window which is exposed to fire products and the insulated or protected edge of the window. 
These stresses lead to cracks and eventually to breakage. In this paper, a thermal analysis of the effect 
of radiation and convection on the glass is carried out by including the exponential decay of thermal 
radiation within the glass and heating and/or cooling of the surfaces due to convection and radiation 


heat transfer. 


Emmons (1986) indicated that the process of the initiation and propagation of cracks in glass due 
to the heat is not very well understood. Finnie et al (1985) carried out an analysis which implied that 
the only factor governing the breakage is the net temperature difference between the central heated por- 
tion of the glass and the protected edge, which can be calculated from the glass properties. Since the 
protected edge temperature is not expected to rise much in a fire, the difference between the tempera- 
ture of the central portion and the edge can be approximated by the net temperature rise of the heated 
portion. The side exposed to the fire is expected have the highest temperature and thus a knowledge of 


the temperature history of only this side is adequate in determining the time for breakage. 


The problem governing this phenomenon is transient and non-linear due to the radiation boundary 
condition. Keski-Rahkonen (1988) linearized the boundary conditions and obtained an exact solution 
for the temperature field. Recently, Davies (1985) used an integral method for determining the tem- 
perature field in a plate losing heat by combined radiation and convection. In Davies’ method, inaccu- 
racies arise due the assumed nature of the temperature profile. For solving the complete non-linear 
problems the use of a numerical method is essential. Methods such as finite-difference, finite-element 


and boundary element have been used before. However the implementation of these methods leads to 


the inaccuracies arising from the discretization of both the time and space domain. Also the whole 
temperature field has to be evaluated. In the present case, the temperature of interest is only of the 


side exposed to the fire and hence the knowledge of the whole temperature field is irrelevant. 


The thickness of the window is usually an order of magnitude smaller than the other dimensions 
and so in the case discussed, the window has been assumed to be an infinite slab. Thus the equation 
governing the system is a one dimensional inhomogeneous heat equation with non-linear boundary 
conditions. To solve this equation, the method of Laplace transform is utilized. After simplifications, 
only two equations for the temperature of each surface need to be solved. This method is similar to 
Chambres (1959) method which was used to determine the temperature distribution in a semi-infinite 
solid with radiation type boundary conditions. In Chambres methods one non-linear Volterra equation 
of the second kind for the surface temperature needs to be solved, whereas in the present case, a sys- 
tem of two coupled non-linear Volterra equations of the second kind for the temperatures of each side 
have to be solved. These equations are solved numerically by using the trapezoidal rule for numerical 
integration and Newton-Raphson’s method for determining the roots of the non-linear equations. The 
equations are exact and inaccuracies arise only due to the discretization of the time domain. This 


method is applied to the general case of time varying radiative flux and surrounding temperatures. 


THEORY AND FORMULATION 


Consider an infinite slab of thickness L, initially at a temperature 7;. At time ¢ = 0 it starts get- 
ting heated on the side x = O due to radiative flux /,/’(¢) and loses or gains heat by combined convec- 
tion and radiation on both sides. Suppose the heat transfer coefficients on the two sides are 
represented by h, and A>, the ambient temperatures are represented by T)..(¢’) and T>..(’ ). Let B be 
the absorption length and & be the thermal concuctivity of glass. Also let p and C,, be the density and 
specific heat respectively of the glass. Also let € and €,, be the emissivities of the glass and the sur- 


roundings. o is the Stefan-Boltzmann constant. Then the governing equation for this system will be 


2 -x'/B 
oT au ar) = (1) 
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where j’ (f) is some specified function of time. The initial and boundary conditions for this problem 


are 


at 2 2 0nT=eh; 


stig 20h 1st = h(T (tT (0, )) - 67’ 4(0,") 
ee oT a , 4 4 
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Nondimensionalization 


This equation is nondimensionalized by using the following dimensionless variables 


with initial and boundary conditions 


Aiath=.0) tia) 
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Ox 
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where 
grt) =A + BT(0,t) + CT2(0,t) + DT7(0,t) + ET*(0,r) 
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METHOD OF SOLUTION 


We use the method of Laplace transforms to solve this equation. The Laplace transform is taken 
with respect to time and is defined as 


oo 


T*= | Te?! dt (20) 


Substituting the transformed variable into the governing equation and applying the boundary conditions 


leads to the solution 
7 elt wae — 
1 Ae apA. ara + B,cosh(Vp x) + C ,cosh(Vp (1-x )) (21) 


where 


A,= (22) 
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We need to solve for both the surface temperatures as the temperature on the hotter surface depends on 
the temperature on the other surface and vice versa. Hence, the transformed form of the temperatures 


are 


T’ (0) = 


Anis = Raya rae as 
: +B,+C,coshYp =f, 4, +f2 49 pet (25) 


—\/y .* 
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Using convolution theorem and after simplifications, the equations to be solved look like 


t t 


(Ou) = fareof evade + [arto fae—ode + [ios se—oat 27) 
(0) 
and 
T (lt) = fame (t—1)dt + [a (tg 9(t—1)dt + 3 [i@este-at (28) 


where f 1, f2, £3, 21, 22, 83 are the kernels which depend upon y and are specified in Appendix I. 
The equations to be solved are coupled non-linear Volterra equations of the second kind. The equa- 


tions are non-linear because both q, and q2 are non-linear functions of 7, and 7 respectively. 


Transformation of variable 
We observe that the integrands are not bounded as f(t) and g,(t) at ¢t =O go to infinity as 
t-'? So if we now transform the equations using the variable u = Vr—t , the equations become 


Nt Vi NI 


T (0,t) = 2 | uF )q\(t-u*)du + 2 | Fol )qx(t—u*)du + - J uF 3(u)j(t—u)du (29) 


and 


Vi Vr VE 
T(1,t) = 2 | uG (ua y(t") du + 2 [uG ru )qx(t—u*)du + 5 [uGaw i (wu (30) 
0) 0 


where F ,(u) = f ,(u) etc. and the integrands are bounded. 


Numerical procedure 


The numerical procedure chosen was trapezoidal rule with constant time steps (thus variable Au ) 
for numerical integration and Newton-Raphson’s methods for finding roots of the non-linear equations. 
By looking at the equations, it appears that the equations are coupled. However, since f ,(0) and g (0) 


are equal to 0, at a particular time step the equations can be solved independently of one another. 


RESULTS 


Figure 1 shows the comparison with the exact solution of the linear problem (given in Appendix 
II) for the case of no heat loss due to thermal radiation. This figure shows the dimensionless tempera- 
ture variation with respect to dimensionless time of the sides exposed and unexposed to incoming radi- 
ative flux. The agreement is within 0.5% for a dimensionless time step of 0.002 suggesting that the 


the time step was quite accurate. 


Figures 2(a) and (b) show respectively the variation with respect to time of the dimensionless 
temperature of the side exposed to and the side unexposed to the incoming heat flux. Only the Biot 
numbers were changed. In this case the initial temperature of the material is set equal to the ambient 
temperatures on either sides. Also the incoming radiative flux is assumed to have a constant value. 
The ambient temperatures are assumed to be constant. 

As expected, the temperature of the side exposed to the heat flux increases more rapidly com- 
pared to the temperature of the unexposed side. As Biot numbers increase, heat loss is greater and 
hence the rate of increase of temperature is lower. 

Figure 3 shows the variation of the temperatures for varying dimensionless absorption length for 


constant value of Biot number. Here also, the input radiative flux function is assumed to be a 


constant. Here, as y decreases, the temperature of the exposed side increases as small y corresponds to 
small absorption length, implying that most of the radiation gets absorbed within a short distance. And 
for large values of y, the material behaves more like a transparent medium. For large values of y an 
artificial temperature drop is encountered as the approximation of surface emission but absorption 


within the body breaks down. 


Figure 4 shows the variation of the temperatures of the sides corresponding to varying input radi- 
ation flux, holding the other parameters constant. It can be seen here that the rate of increase of the 
temperature of the exposed side varies roughly as the integral with respect to time of the incoming 
radiation function. This is because the rate of increase of temperature is dependent on the rate of 


incoming heat as can be seen in equation (29)for T (0,r) . 


Table 1 shows the various values used for the calculations. These values represent the typical 
values which might be encountered in the situations of real fires. The time of glass breakage solely 
depends upon the net temperature rise of the glass, and it was found that for ordinary plate glass, the 
value was found to be about 60° C. This value is obtained from Hooke’s law and it was verified by 
Finnie et. al. that this value does not depend on the rate of temperature rise. This would imply that 
for a particular value of the radiative flux, a value of the dimensionless critical temperature can be 
obtained. Hence in the computer program, the computations are stopped as soon as the cnitical tem- 


perature is reached. 
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APPENDIX I 


The Laplace transforms of the kemels in equations are 
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APPENDIX II 


The exact solution to the linear problem of negligible heat loss due to radiation is 


T=u(%,t) + v(x) (45) 
where 
v(x) = px +q (46) 
where 
ABi, + BBi, B -—A —ABi, 
a ae (47) 
Bi, + Bi, +Bi, Bi, Bi, + Bi +Bi Bi, 
Here 
pd Gs ae eA es oy 
A= Biba (48) 
Te ; 1P 
where 
h,L hob 
Bi, = — and Bi» = (49) 
k k 
where h, and hp are the heat transfer coefficients. The solution for u(x,t) is given by 
co no) t ee | 
We Syie sees fe™ (Oy (t)dt| 0, (x) (50) 
n=0 
4 
where the characteristic functions are 
o, (x) =A, cos(A, x) + Bizsin(A, x) (S1) 
and the eigenvalues A, are obtained from the solutions to 
hee BitBi 
COL Neal eee eae (52) 


2, (Bi }+Bi>) 


The constants C, appearing in the equation for u(x,t) are defined as 


1 
|»), dx (53) 


1 
N(A,) 


the weight functions w, are defined as 


1 


Ww, (t) = Lime *%p, (x)dx 


1 
YN (A, ) 


and norms N (A, ) are defined as 


N(A,) = 


No 


he 


(A,° + Biz?) 1 + —=——| + Bi, 
A, + Bi, 


(54) 


(55) 


thermal conductivity, 


thermal diffusivity, 


thickness, 


penetration depth, 


temperatures, 


emissivities, 


incident flux, 


Table 1. Glass and Fire Parameters 


k = 0.76W /mK 

= 3.60 10m nts 

L = 6.25 X 107m 

B=10%m, 3.13X10°m, 6.25X 10% m 


T 0 = 12. = T; = 293 K 


I) = 1000 W/m?K 
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FIGURE CAPTIONS 


Figure 1. The surface temperatures, T7(0,t) and 7(1,t) are shown where T = (T’-T;)k/I9£ and 
t =r o/L* for y= 0.16 and Bi, = Bi, = 0.1. The exact linearized results and numerical results are 


indistinguishable. 


Figure 2. The variation of the exposed surface temperature with Biot number for Bi = Bi, = Biz is 


shown for y = 0.16 and the other parameters as listed in Table 1. 


Figure 2. The variation of the unexposed surface temperature with Biot number for Bi = Bi, = Bi, is 


shown for y = 0.16 and the other parameters as listed in Table 1. 


Figure 3. The variation of the exposed surface temperature with depth parameter, y, for Bi, = Biz = 0.1 


and other parameters listed in Table 1. 


Figure 4. The effect of different time dependencies of the incident radiative flux on the exposed surface 


temperature for y = 0.16, Bi, =Bi. = 0.1 and the other parameters as listed in Table 1. 
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